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Abstract 

We measured the single-photon detection efficiency of NbN superconducting single photon de- 
tectors as a function of the polarization state of the incident light for different wavelengths in the 
range from 488 nm to 1550 nm. The polarization contrast varies from ~5% at 488 nm to ~30% at 
1550 nm, in good agreement with numerical calculations. We use an optical-impedance model to 
describe the absorption for polarization parallel to the wires of the detector. For lossy NbN films, 
the absorption can be kept constant by keeping the product of layer thickness and filling factor 
constant. As a consequence, we find that the maximum possible absorption is independent of filling 
factor. By illuminating the detector through the substrate, an absorption efficiency of ~ 70% can 
be reached for a detector on Si or GaAs, without the need for an optical cavity. 



* Electronic address: 



driessen@molphys.leidenuniv.nl 



I. INTRODUCTION 



Superconducting single-photon detectors (SSPDs) [T], that consist of a meandering NbN 
wire, are an interesting new class of detectors that may outperform single-photon counting 
avalanche photodiodes. SSPDs feature a relatively high quantum efficiency at infrared wave- 
lengths, combined with low time jitter, low dark counts, and high counting rates [2|. This 
makes these detectors promising for quantum optical studies and long-distance quantum 
cryptography applications [3]. 

A lot of attention has been given to the electronic operation of these detectors [H 0, E], 
leaving the optical design of the detectors less explored. In fact, due to the highly anisotropic 
nature of the wires, the detection efficiency shows a strong polarization dependence [7]. This 
is important, since a common way to encode quantum information is to use the polarization 
state of the photons |8]. Detection of a photon thus comprises a simultaneous measurement 
of the polarization, which may be undesirable for some apphcations. At the same time, 
knowledge of the polarization dependence may simplify experimental schemes that require 
a polarization measurement, or can be used to optimize the detection efficiency. 

The efficiency rj to detect a single photon can be decomposed in an electronic and an 
optical contribution and can be expressed as 

V = VeA (1) 

where A is the optical absorption efficiency of the detector, and r/e is the electronic efficiency 
of the detector, i.e. the probability that an absorbed photon leads to a measurable voltage 
pulse across the detector. 

The microscopic working principle of the detectors, which is essential to understand r/e, is 
still under active investigation [9l[l0]. On a macroscopic level, a photon that is absorbed by 
the superconducting wire triggers a temporary loss of superconductivity, which gives rise to 
a finite voltage pulse across the detector. The optical absorption efficiency A is determined 
by the geometry of the detector and the dielectric constants of the substrate and the NbN 
layer. Since the energy of the incident photons is much larger than the superconducting gap 
of the NbN, the complex dielectric constant of the NbN layer at room temperature can be 
used. 

The polarization dependence of NbN SSPDs has been investigated at a single wavelength 



and compared to finite-difference time domain calculations f^. In section III, we experi 



mentally investigate the wavelength dependence of the polarization contrast, in the range 
between 488 and 1550 nm, and report a strong dependence of the polarization contrast on 



the wavelength. We introduce an analytical optical impedance model in section [IV] to de- 
scribe the optical absorption in thin lossy films and describe different ways to increase the 
detection efficiency by changing the parameters of the detector. We find that the optimum 
thickness is a strong function of the fill fraction, while the maximum achievable absorption 
is independent of the NbN fill fraction. 

The optical impedance model also provides more insight into the cavity enhancement re- 



ported for a NbN detector inside a Fabry-Perot type cavity [3 HI]. We show, in section IV C 



that the absorption of the detector is enhanced by a factor n, with n the refractive index of 
the substrate, when the detector is illuminated from the substrate. This factor was not ac- 
counted for in earlier work and thus leads to an overestimate of the resonant enhancement. 
For a high index Si or GaAs substrate this factor becomes dominant and an absorption 
efficiencie of ~ 70% can be reached without the need of an optical cavity. 



II. EXPERIMENTAL SETUP 

In our experiments, we used a commercial NbN SSPD [2], with an area of 10x10 /im^. 
The detector consists of a ~4 nm thick NbN meander on a R-plane sapphire substrate. It 
has a nominal fine width of 100 nm and a filling factor of ~ 55%. Fig. [T](a) shows a scanning 
electron microscope (SEM) image of a detector similar to the one used in the measurements. 

We mounted the SSPD in a ^He-cryostat and cooled it to a temperature of ~5 K. The 
temperature remained constant within 10 mK during each measurement run. Figure [ij^b) 
shows a schematic overview of the electronic circuit used to operate the detector. The 
detector was biased at 90% of the critical current through a bias-T with a 400 kfi resistor. 
The equivalent circuit of the detector (dashed box) contains a switch that is closed in the 
superconducting state. When a photon is absorbed, the switch opens temporarily [5]. The 
resulting voltage pulse across the detector is amplified (66 dB) and detected by pulse counting 
electronics. 

Unpolarized light from an incandescent tungsten lamp was wavelength-filtered and sent 
through a 50 /im core size multimode optical fiber. The output of the fiber was imaged onto 
the detector using a telescope and a lens mounted on a piezo stage inside the cryostat, as 
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FIG. 1: (a) SEM image of a NbN SSPD similar to the one used in the experiments. The 100 nm 
wide NbN line is folded into a meander with an area of 10x10 ^m^. (b) Schematic diagram of 
the readout electronics. The bias current is provided by a voltage source and a resistor of 400 kO. 
The dashed box contains a phenomenologically equivalent circuit of the detector, (c) Schematic 
overview of the optical setup. Wavelength-filtered light from a lamp is sent through an optical 
fiber with a 50 fim core, and is imaged onto the detector through a telescope and a moveable lens 
mounted inside the cryostat. 

shown in Fig. [T|^c). To probe the polarization dependence of the detection efficiency, a linear 
polarizer with an extinction ratio better than 100 : 1 for the wavelength range of interest 
was placed in the parallel part of the beam. To probe the wavelength dependence, we used 
different narrow bandpass filters (<10 nm FWHM) in combination with several edge filters 
to ensure that the fight on the detector was monochromatic. 

III. POLARIZATION DEPENDENCE 

Figure [2] shows the count rate of the detector as a function of linear polarization for a 
wavelength of 1550 nm (black squares) and 532 nm (red triangles). Note that the absolute 
count rates at different wavelengths cannot be compared directly, due to a difference in 
incident power. The insets show the orientation of the £'-field relative to the detector. The 
measured count rates follow a sinusoidal dependence as a function of polarization and are 
minimal when the E-&eld is perpendicular to the lines of the detector. 
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We define the polarization contrast C as 



C 



iV|| - Nx 
iV|| + 



(2) 



where A''|j and A^^ are the count rates of the detector when the light is polarized parallel and 
perpendicular to the wires, respectively. This definition of the contrast is a direct measure 
for the polarization effects, independent of the electronic quantum efficiency (77e), and the 
incident power. We extract the contrast from the sinusoidal fits to the data (solid curves in 
Fig. [2]). It varies with the wavelength of the incident light and is independent of the bias 
current and temperature of the detector in our experiment. 
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FIG. 2: Count rate of the SSPD (corrected for stray light counts) as a function of linear polarization, 
for a wavelength of 1550 nm (black squares) and 532 nm (red triangles). The blue circles show the 
count rates for 1550 nm light when two wedge depolarizers (under a relative angle of 45°) are placed 
after the polarizer. The insets show the orientation of the i?-field relative to the detector, for the 
different polarizer settings. The soHd curves are sinusoidal fits to the data, used to extract the 
polarization contrast. 

The blue circles in Fig. [2] show the count rate as a function of polarizer angle, at a 
wavelength of 1550 nm, when two wedge depolarizers under a relative angle of 45° were 
placed after the polarizer. These wedge depolarizers effectively depolarize the incident light 
by imposing a position-dependent rotation of the polarization. Indeed, the polarization 
contrast in this case is reduced to below 3%. The lower average count rate can be attributed 
to the extra four air-glass interfaces in the optical setup, leading to an increased refiection 
of the incident light. 



The polarization effect can be understood by comparing the periodic structure of the 
detector to that of a wire grid polarizer |T2] that consists of a grid of parallel, highly conduc- 
tive metal wires with a subwavelength spacing. For a perfect conductor the E-iie\d should 
be perpendicular to the metal surface. As a consequence, only light with a polarization 
perpendicular to the wires is efficiently transmitted. A similar argument holds for lossy 
metals, albeit that in this case the field penetrates into the metal, leading to absorption. 
This absorption is largest when the E-Held is parallel to the wires, since in this case the field 
penetrates more into the metal. 

For the typical dimensions and spacing of the NbN wires, an effective medium approach 
that is accurate for both polarizations is difficult [13], [H]. Instead, we calculated the ab- 
sorption at normal incidence for an infinitely-sized detector, using the rigorous coupled-wave 
analysis (RCWA) developed in Ref. 14. This method finds an exact solution of Maxwell's 
equations by expressing the electromagnetic fields in the different materials as a summa- 
tion over all diffraction orders. The Fourier components of the periodic dielectric constant 
couple the diffraction orders in the patterned region. The continuity of the parallel compo- 
nent of the wavevector, together with the boundary conditions for the E and H fields fully 
determine the field in all regions. From this the intensity in all reflected and transmitted 
diffraction orders can be calculated. The absorption in the grating is then simply given by 
A = 1 — R — T, where R and T are the reflected and transmitted intensity. 

The effects of focusing of the incident beam can be taken into account by decomposing 
the beam into plane waves with wave vector k. Each of these plane waves will experience a 
different absorption A{k). The effect of finite detector size can be incorporated in a similar 
way, by multiplying the beam profile in the near field by an aperture function D{f) which 
is 1 at the location of the detector, and elsewhere. Taking both into account, the total 
absorption is given by the convolution integral 



of the aperture function D{r). 

The /c-spread of the incident waves is determined by the detector size [determining the 




(3) 



where u{k) is the Fourier transform of the beam profile, and D{k) is the Fourier transform 



spread in D{k)] and the numerical aperture of the last lens in the illuminating system, 
determining the spread in u{k). The latter is the most important factor in our experiment. 
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since we used a large-NA lens to focus the incoming light onto the detector. Calculations 
of the absorption of the grating as a function of angle of incidence ( 1.6., ctS ct function of k) 
show however, that the absorption only varies appreciably from the absorption at normal 
incidence for angles of incidence corresponding to NA > 0.5. Therefore, the total absorption 
given by Eq. ^ can be approximated by a product of the absorption coefficient at normal 
incidence and the total intensity impinging on the (finite-sized) detector. This justifies the 
use of a plane wave calculation in the rest of this Paper. 

To calculate the absorption efficiency, we used the nominal structure parameters of the 
detector, and tabulated values of the dielectric constant of the sapphire substrate (^sapphire = 
1.74 at 1550 nm) [16]. For the wavelength-dependent dielectric constant of NbN, a Drude 
model [IT] was used, giving a refractive index nNbN = 5.5 + 6. 3z at a wavelength of 1550 nm. 
This value is close to the value reported in Ref. 7, for a thicker NbN film. 

Figure [3|a) shows the calculated absorption for polarization parallel (blue fine) and per- 
pendicular (red line) to the wires, as a function of wavelength. The absorption for parallel- 
polarized light monotonously increases with wavelength, whereas the absorption for perpen- 
dicular polarization goes through a maximum and decreases for wavelengths above 800 nm. 
This leads to a higher polarization contrast for longer wavelengths. 

For comparison, the dashed line in Fig. [3|a) shows the absorption of an unpatterned film, 
multiplied by the filling factor of NbN, as was suggested in Ref. 1. This estimate deviates 
over the entire wavelength range from the polarization-averaged result obtained by RCWA, 
which shows that for structures with features smaller than the wavelength of light, a more 
refined model is needed. We will discuss this refined model in Sec. HVl The fact that the 
absorption decreases for both the parallel polarization and for the closed film is mostly due 
to dispersion of the dielectric constant of the NbN material, eNbN- 

In Fig.|3|b) we compare the measured polarization contrast (red dots) to the results of the 
calculations (black solid curve), as a function of wavelength. For comparison, the calculated 
contrast is shown for filhng factors of 52% (dashed curve) and 58% (dash-dotted curve) as 
well. The experimentally observed contrast varies between ~5% and ~30% and increases 
with wavelength. The error bars on the experimental points represent slight variations in 
the measured polarization contrast during different measurement runs, as well as a slight 
polarization in the illuminating fight source, of ~ 1%. We attribute the fact that the 
calculation and the measurements differ for lower wavelengths to the fact that we used 
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FIG. 3: (a) Calculated absorption efficiency of a NbN grating as function of wavelength for po- 
larization parallel (blue curve) and perpendicular (red curve) to the lines of the detector. For 
comparison, the dashed line shows the calculated absorption of an unpatterned film multiplied by 
the fining factor of NbN. (b) Measured (red dots) and calculated (black curves) polarization con- 
trast as a function of wavelength. The calculations are shown for a fiUing factor of 52% (dashed), 
55% (solid), and 58% (dash-dotted) and a film thickness of 4.5 nm. 

literature values for the dielectric constant of NbN. It is known that the dielectric constant 
of NbN varies as a function of the deposition parameters [17] and may depend on the film 
thickness as well [Hj. Additional calculations (not shown) reveal, that for lower wavelengths, 
the polarization contrast is increasingly sensitive to small variations in the dielectric constant 
of NbN. 

It has been shown that the linear-polarization dependence can be removed by changing 
the design of the detector [19]. A spirafing detector breaks the translational symmetry that 
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causes the strong polarization contrast. The optical absorption in such a detector, however, 
will be lower than the maximum obtainable for parallel-polarized light, due to the fact that 
in these detectors, partial screening of the electric field is always possible. 



IV. DISCUSSION 

A. An optical impedance model for the absorption of a metal film 

In order to gain some physical insight into the absorption in the detector, we start out 
by describing the absorption of a film of thickness d with a complex dielectric constant €2, 
embedded between two dielectrics with refractive index rii and ns, respectively. The film is 
illuminated from the medium with index ni. 

We can define the optical impedance of a medium i with refractive index rii as 

m = ^. (4) 

rii 

where rjo = \/ijo/^o = 377n is the impedance of the vacuum. The reflection and transmission 
of the layered system are given by [20] 

2 



R 



Vload - Vl 



T = ^ 

^3 



^?ioad + m 



(5) 
(6) 



'7load + Vl 

where T^ioad is the combined load impedance of the film and the backing substrate. The 
absorption of the film is again given hyA = l — R — T. 

If we assume that the film is thin enough to neglect interference effects (kod <^ 1), the 
load impedance is given by [2lj 

Rn + V3 

where Rn ~ rjo/kod Im 62 is the square resistance for a highly absorbing (Im €2 ^ Re €2) 
film, and ko is the wave vector of the light in vacuo. With these assumptions, we can write 
the absorption in the film as 

A=^(^I1^^^X = 4n, ^-^^^^ (8) 

ViRa \Vi + Ra + J {ni + + kod Im €2) 
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The absorption of the film reaches a maximum value Amax = ni/{ni + n^) for a square 
resistance given by 

Rn = (9) 

rii + ?73 

Note that the maximum possible absorption is a function of the refractive indices of the 
surrounding media only. The optimal value of Ru to reach this maximum can be obtained 
by tuning the film thickness d. 

B. The effect of film thickness 

Figure [4] shows the absorption and the polarization contrast of a film of NbN, embedded 
between air (rii = 1) and sapphire (^3 = 1.74), as a function of the film thickness. The solid 
curves show the calculated absorption using the rigorous coupled-wave analysis described 
before, while the dotted curves are obtained from the impedance model. 

For a closed film (black curves), there is a distinct maximum of absorption, that occurs 
at a thickness 

ko Im 62 ' 

For thinner films, the transmission through the film is too high to get maximal absorption, 
whereas for thicker films, refiection dominates. 

The blue and red curves in Fig. [4] show the absorption for a detector with filling factor 
0.5 and lattice period 200 nm, for polarization parallel and perpendicular to the wires, 
respectively. The thickness for which the absorption in the patterned film is maximum, is 
higher than the optimal thickness for the closed film. The dotted line is calculated using the 
impedance model of section |IV A| taking an effective dielectric constant for the absorbing 



film, given by [13] 

Ccff = (1 ^ /)eslits + /cNbN, (11) 

where / is the filling factor of the metal, and esuts is the dielectric constant of the material 
in the slits, typically air. Since only the imaginary part of ecfr determines the absorption 
in the film, the absorption of the detector can simply be calculated by multiplying the 
thickness of the film by the filling factor. For the polarization perpendicular to the wires 
of the detector, it is not so straightforward to define an effective dielectric constant for the 
patterned film [131 [H]. For this polarization the fight is concentrated in the air slits and the 



10 




Thickness (nm) 

FIG. 4: Calculated absorption at a wavelength of 1550 nm as a function of film thickness. The 
black curve gives the absorption for a closed film, the blue and red curves for a detector with lattice 
period 200 nm and filling factor 0.5, for polarization parallel (blue) and perpendicular (red) to the 
wires of the detector. The dotted curves are calculated using the impedance model of section [TVX 
whereas the solid curves are exact calculations using RCWA. The top graph shows the calculated 
polarization contrast. 

effective dielectric constant is closer to that of air. Therefore the condition Im 62 3> Re 62 
used to define the impedance model, breaks down for this polarization. 

Surprisingly, the calculation also shows that it is easily possible to construct a detector 
where the absorption for parallel polarization is larger than the absorption of an unpatterned 
film of the same thickness. Since the electronic efficiency of the detector, T^e, strongly depends 
on the thickness of the metal [22l [23], it is important to realize that the absorption for 
parallel-polarized light is a function of df Im 62. A reduction in thickness of the detector, 
to increase the electronic efficiency, can thus be countered by increasing the filling factor 
accordingly. 

C. Illuminating through sub- or superstrate 

Commonly, NbN SSPDs are deposited on a substrate of sapphire and illuminated from 
air. An inspection of Eq. ([s]) shows that for a certain choice of sub- and superstrate, a factor 
of n^/rii in absorption can be gained by illuminating the detector from the medium with 
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the higher refractive index. Figure [5] shows the calculated absorption for a detector, with 
a superstrate of air (ni = 1), as a function of the refractive index of the substrate. The 
thickness of the detector is set such that maximal absorption in the detector is achieved. 
This thickness is indicated with the black line. The solid curves give the absorption for 
illumination from the air, whereas the dash-dotted curves give the absorption for illumination 
from the substrate. The blue and red curves are for polarization parallel and perpendicular 
to the wires, respectively. We stress that this effect is caused by a lower impedance mismatch 
and should be separated from the cavity enhancement of the absorption, previously reported 
in Refs. [3 [II]. 




Substrate index 



FIG. 5: Calculated absorption for a detector with filling factor 0.5, lattice period 200 nm, at a 
wavelength of 1550 nm, as a function of the substrate refractive index. The solid curves are for 
illumination from the air side, the dash-dotted curves for illumination from the substrate side. The 
blue and red curves give the absorption for polarization parallel and perpendicular to the wires, 
respectively. The detector thickness is changed at each substrate index, to achieve maximal absorp- 
tion. The thickness is given by the black line (right axis). The top graph shows the polarization 
contrast. 

When the substrate index is increased, the absorption rises for iUumination from the 
substrate side. For illumination from the air side, the absorption for parallel polarization 
decreases. Note however that in both cases, the polarization contrast decreases, from C = 
0.88 at na = 1 to C = 0.26 at ns = 4, and is independent on the direction of iUumination, as 
shown in the top graph of Fig. [5| The absorption is a factor of higher, when the detector 
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is illuminated from the substrate, as expected from the impedance model. It is interesting to 
note that, for parallel-polarized hght, the absorptions from super- and substrate add up to 
give Auper+^sub ~ 1- It is therefore possible to construct a detector with higher absorption, 
up to 70%, and lower polarization contrast, by using a high refractive index substrate (e.g. 
Si or GaAs) and illuminating the detector from the substrate. Unfortunately, increasing the 
refractive index of the substrate also increases the wavelength for which diffraction orders in 
the substrate appear. The first diffraction order at normal incidence appears at X/a — ris, 
with a the periodicity of the structure, and A the wavelength of the light. In general, these 
diffraction orders lower the absorption efficiency. For a typical lattice period of 200 nm, and 
a substrate index of = 3.5, the first diffraction order appears at a wavelength of 700 nm, 
making detectors on a high-refractivc-indcx substrate less effective for detecting visible light. 
The problem of diffraction could also be circumvented by designing a detector that has a 
variable line spacing. 



V. CONCLUSION 



In conclusion, we have measured a polarization dependence in the detection efficiency of 
NbN superconducting single photon detectors and find a wavelength dependent polarization 
contrast between 5% and 30%. This effect can be explained by the geometry of the detector. 
Calculations of the optical absorption efficiency give good agreement with the measured data. 
We have demonstrated that the polarization dependence can be removed by the use of wedge 
depolarizers. 

Furthermore, we have shown that the parameters of the detector can be tuned to achieve 
an absorption for a polarization parallel to the detector wires, that exceeds the absorption 
of an unpatterned film of the same thickness. We have given a simple optical impedance 
model, that allows for a quick estimate of the parameters needed to optimize the detector. 
For parallel-polarized light, the maximum absorption achievable is not determined by the 
thickness or the dielectric constant of the metal film, nor by the filling factor, but only by 
the refractive indices of the surrounding media. We have shown that by illuminating the 
detector from the substrate it is possible to increase the detection efficiency of the detector 
even further, by a factor equal to the refractive index of the substrate. Such highly absorbing, 
highly polarization-dependent detectors can be employed to efficiently detect photons with 
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a well-defined polarization. 
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